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The Hong-Ou-Mandel (HOM) effect is widely regarded as the quintessential quantum interference 
phenomenon in optics. In this work we examine how nonlinearity can smear statistical photon 
bunching in the HOM interferometer. We model both the nonlinearity and a balanced beam splitter 
with a single two-level system and calculate a finite probability of anti-bunching arising in this 
geometry. We thus argue that the presence of such nonlinearity would reduce the visibility in the 
standard HOM setup, offering some explanation for the diminution of the HOM visibility observed 
in many experiments. We use the same model to show that the nonlinearity affects a resonant two- 
photon propagation through a two-level impurity in a waveguide due to a “weak photon blockade” 
caused by the impossibility of double-occupancy and argue that this effect might be stronger for 
multi-photon propagation. 

PACS numbers: 42.50.-p, 42.50.Ar, 42.50.Ct 


INTRODUCTION 

The Hong-Ou-Mandel (HOM) interferometer pQ [2] is 
one of the main tools in registering biphotons (i.e. entan¬ 
gled photon pairs) created by SPDC (spontaneous para¬ 
metric down-conversion ) pHTTl . The ideal HOM protocol 
comprises a four-port interferometer with two incoming 
and two outgoing channels, a spectrally-flat, balanced 
beam splitter and the coincidence counter (see Fig. [T]) . 
When uncorrelated photons arrive at the incoming port, 
the outgoing photons split equally between the two detec¬ 
tors resulting in a signal in the coincidence counter. How¬ 
ever, when two photons arrive simultaneously at both 
the incoming ports, the outgoing photons are bunched 
together leaving the interferometer through only one of 
the outgoing ports mm, so that the counter detects no 
signal. As the entangled photons are simultaneously cre¬ 
ated by SPDC, their arrival at the two incoming ports 
(against the background of uncorrelated arrivals of non- 
entangled photons) is manifested by a dip in a signal 
registered by the coincidence counter. 

In this work we argue that characteristics of such a dip 
can be substantially affected by the presence of a nonlin¬ 
earity that leads to an effective interaction between two 
photons that simultaneously enter the beam splitter in 
the HOM device. The nonlinearity might be due to a 
single two-level impurity embedded in the device. The 
interplay between such a nonlinearity and HOM interfer¬ 
ence leads to a novel mechanism of the suppression of the 
photon bunching which might obscure information about 
biphoton generation. 

An alternative geometry where nonlinearity is also es¬ 
sential is illustrated in Fig. [U*- a two-level impurity is 
embedded in a waveguide. Such an impurity would prac¬ 
tically not affect the propagation of off-resonance photons 
but result in a full reflection of a single photon at the 
resonant frequency, acting similarly to the side-attached 
resonant impurity for electrons propagating through a 



FIG. 1. Traditional Hong-Ou-Mandel interference scheme. 
Two identical photons arriving simultaneously at a balanced, 
broadband beam splitter (BS) will be conveyed along only 
one of the possible outgoing channels - and so, in contrast to 
the general case, no coincidence counts will be accumulated. 

conducting wire j!2j . In this geometry the nonlinear¬ 
ity suppresses the resonant reflection and leads to a par¬ 
tial propagation for the photon pair at the resonance. A 
feasible “dual” geometry, Fig. ih comprises two waveg¬ 
uides separated by an optically opaque region with an 
embedded two-level impurity. There is no propagation 
of non-resonant photons through such a weak link but 
an ideal resonant transmission of a single photon, which 
is again suppressed by the nonlinearity. We show that all 
such geometries can be described within essentially the 
same model where the nonlinearity substantially alters 
two-photon propagation. 

THE MODEL 

We consider a wave-guided few-photon beam interact¬ 
ing with a single near-resonant atom that can be de¬ 
scribed as a two-level system (TLS). We assume that 
both incident and transmitted or reflected photons can 
propagate along two channels. Such a model mimics (un- 
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FIG. 2. Two additional geometries described by the model 
of Eq. 0: (a) a single atom (TLS) embedded in a ID pho¬ 
tonic crystal waveguide leads to reflection of resonant photons 
in the channel (off-resonant photons are freely transmitted); 
(b) an interstitial TLS provides a resonant link between two 
waveguide channels. 


der conditions specified below) the HOM geometry of 
Fig. [TJwhen the incoming channels are different, and cor¬ 
responds to a TLS embedded into a ID photonic crystal 
waveguide when photons enter through the same chan¬ 
nel, Fig.[2j The corresponding Hamiltonian is reduced in 
the usual rotating wave approximation to 
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Here b+ k and b a k are the photon creation and annihi¬ 
lation operators with the index a = 1,2 labelling the 
channels; ui = k is the photon energy (counted from the 
upper energy level of the TLS) in the units where both 
h and the group velocity of light in the medium are set 
to 1; the TLS is described by Pauli’s rasing and lowering 
matrices, a± = a x ± ia y ; and y/j is the atom-photon 
coupling strength, with 7 being the TLS relaxation rate. 
Neglecting the nonlinearity, i.e. considering two- or multi¬ 
photon scattering from the TLS as totally independent, 
this model describes - with a proper choice of 7 and w 
specified below - both a 50:50 (balanced) beam splitter 
for the HOM geometry and ideal resonant reflection or 
transmission for the geometry of Fig. [2ja) and (b), re¬ 
spectively. 

We will describe how the effects of the nonlinearity 
caused by indirect photon interaction via scattering from 
the TLS change the probability P aa '\ppi of two photons 
entering via channels /3 and /3' and exiting via channels 


a and a' in the Results below. In both cases of the HOM 
geometry (/3 ^ /?') and the resonance geometry (/3 = /?') 
we calculate the probability of photons leaving through 
different channels: Prom = Pi 2 |i 2 + P 21112 and P res = 
P\ 2 \pp + P 2 i\pp- The former describes the suppression of 
photon bunching in the HOM experiment while the latter 
describes a “weak photon blockade” of the resonant two- 
photon reflection from or transmission through TLS. 

Before describing the nonlinearity effects, we express 
the probability -P^a' 1 / 9 / 3 ' °f independent two-photon scat¬ 
tering from the TLS via parameters of the model. One- 
photon scattering from a TLS is described by the scatter¬ 
ing matrix S = {s a / 3 } connecting incident and outgoing 
photons in the two channels: 

Hout = S a/3 I0>in > S=(^), ^ = i- (2) 


The unitarity of S-matrix requires the equality of the 
reflection and transmission amplitudes, \t\ = \t’\ and 
|r| = |r'|. In a system with reflection symmetry the 
choice r = r' and t = t' is assumed in what follows. 

For a monochromatic photon in the resonant transmis¬ 
sion geometry of Fig. p'b), the transmission and reflec¬ 
tion amplitudes in Eq. (21) become 


LJ 
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where ui is a photon detuning from the resonance, and the 
TLS is assumed to be initially in its ground state. The 
dual geometry of Fig. [2ja) is obtained by swapping the 
reflection and transmission amplitudes, r ft J. In both 
these cases (choosing the incoming channel /3 = j3' = 1), 
P° u is just a product of one-photon probabilities so 
that Py t |= 1 in the geometry of Fig. |2|a) (both the 
photons are resonantly reflected), or P° 2 in = 1 dr the 
dual geometry of Fig. §b) (both are resonantly trans¬ 
mitted). Only nonlinearity will lead to distinct features 
in two-photon scattering (as compared to single-photon) 
in this geometry. 

For the HOM geometry (/? ^ /T), the probability of 
anti-bunching without nonlinearity is 


Phom — l s n s 22 + S 12 S 21 I 2 — 1 1 2 + r 2 | 2 — (T — P) 2 , (4) 


where T = \t\ 2 and R = |r | 2 are the single-photon 
transmission and reflection probabilities. Thus a perfect 
bunching, Prom = 0, occurs when T = R = \ so that 
the model of Eq. 0 would emulate an ideally balanced 
beam splitter when the frequency ui of monochromatic 
incoming photons coincides with the resonance width 7 . 
For realistic time-resolved photons, an almost balanced 
beam splitter will be emulated if both photons have spec¬ 
tral functions centred at u> = 7 of width a much smaller 
than 7 . The limitation r = 7 /er>lis essential for the 
model we consider here. 
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FIG. 3. The anti-bunching probability in the HOM geom¬ 
etry. Neglecting the nonlinearity (independent scattering), 
it goes to 0 with r —► oo (the monochromatic limit), which 
corresponds to the balanced beam splitter. The nonlinearity 
makes the probability of anti-bunching finite for any value of 
T (the solid line describes the exact numerical solution); the 
analytic asymptotics (the strong-coupling - monochromatic - 
limit) practically coincides with the exact solution for T > 2. 



FIG. 4. Weak photon blockade. Here <5 P leE is the nonlinearity- 
induced anti-bunching probability in the resonance geome¬ 
try of Fig. 0b), i.e. the probability of reflecting one out of 
two photons, as a function of their relative delay A = <jt. 
The inset shows a non-monotonic behaviour of the total anti¬ 
bunching probability due to a trivial off-resonance reflection 
at smaller F. 


This result holds for two identical photons. If the 
photons become distinguishable, e.g. by a delay r be¬ 
tween their arrival times, they would be uncorrelated for 
r a^ 1 and thus have equal probabilities to go to ports 
1 or 2. Therefore, one expects a dip in Phom( t ) at t = 0 
with a non-universal shape dependent on the spectral 
function. When the latter is Gaussian, one has 


Phom( t ) 



(5) 


The depth of such a dip is thus used, e.g., to characterise 
a rate of the SPDC biphoton production, as the shape of 
the entangled photons is identical and they are created 
simultaneously. 

Here we argue that a non-ideal dip may result from the 
nonlinearity in the HOM beam splitter like that described 
by the Hamiltonian 0 - 


RESULTS 

In the HOM geometry, when the identical photons 
come from the opposite (/3 ^ /3') channels without a time 
delay, the nonlinearity results in a nonzero probability 
-Phom(t 7 ^ 0 ) of detecting the two photons in different 
outgoing channels. We find the probability of the coinci¬ 
dence counter clicking as a function of the time delay r 
between two photons hitting the TLS as 

P’hom('t) = Thom( t ) + ^-Phom('t) , (6) 

where the last term corresponds to the nonlinearity- 
induced photon-photon interaction. In particular, for 


photons with narrow Gaussian spectral function |13| . 
when T» 1 , we find 

<5-Phom(V) = —L—e -A , A = err, T = —. (7) 

sJtxY a 

Since ^Thom(O) does not vanish (which is true for any 
photon spectral function), the HOM dip is no longer 
ideal. A more general case, not restricted to T 1, 
is illustrated in Fig. [3] while general expressions for both 
factors in Eq. ([ 6 | will be given in Methods below. 

In the geometries of Fig. [2] two incident photons tuned 
into the resonance w = 0 are coming through the same 
channel (/3 = (}' = 1 ). As the two geometries are dual, 
we give the results on that of Fig. In this case the 
nonlinearity suppresses the resonant transmission as two 
photons cannot simultaneously excite the TLS. Such a 
“weak photon blockade” by TLS results in the nonzero 
probability, P res , of one photon transmitted and the other 
reflected. This probability is peaked for simultaneous 
photon arrival, A = 0. For photons with narrow Gaus¬ 
sian spectral function we find 

<5-Pres = Pres — P r es = Cbl = l = f^ 7T ) (8) 

y 7Tl 1 + e^ 

In this limit the blockade is weak even in the absence of 
the time delay r since the TLS re-emits photon over the 
time interval q ^ 1 which is much shorter than the time 
interval a^ 1 over which each photons arrives. However, 
with r decreasing both the nonlinearity becomes weaker 
and a trivial single-photon off-resonant scattering plays a 
bigger role. Thus the full anti-bunching probability has 
a maximum at T ~ 1 as shown on the inset in Fig. [4j 
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METHODS 


A two-photon incoming state can be represented as 


The multi-photon scattering states for the model 0 
can be described in terms of the Bethe ansatz eigenstates, 
see e.g. Refs. mug, each state parameterised by a string 
of complex quasi-nromenta. The physical implication of 
these complex-valued quasi-momenta is the emergence 
of many-photon states weakly bound to the atom (TLS). 
The asymptotic scattering states constructed out of these 
eigenstates will be affected by such binding. 

For the two-photon scattering states we will use a 
more straightforward and physically transparent ap¬ 
proach based on the scattering matrix. Note that in 
two-photon (or multi-photon) scattering described by the 
model 0 two opposing tendencies compete: the statisti¬ 
cal photon bunching clearly demonstrated by the dip in 
the coincidence counter in the idealised HOM geometry, 
and the nonlinearity caused by the fact that at any mo¬ 
ment only a single photon can go through the TLS. It is 
this competition that spoils the ideal antibunching in the 
HOM geometry and the ideal resonant transmission (or 
reflection) in the geometries of Fig. [2] 

Scattering properties of the atom-light interaction de¬ 
scribed by Hamiltonian (jT]) are encoded in the S'-matrix, 
Eqs. ([2])-(|3|, which links asymptotic incoming states to 
the corresponding outgoing states. Inevitable sponta¬ 
neous emission means that in both the incoming and out¬ 
going asymptotic state the TLS is in its ground state so 
that in specifying the asymptotic states one needs to re¬ 
fer explicitly only to their photonic part. To allow for 
time-resolved photons, in what follows we consider pho¬ 
tons with a spectral function g{oj — wo), centred around 
uiq with a spectral width a. For analytic calculations it 
can be chosen Gaussian, 


= j \P,q-,P',q') (10) 

where Bpp'(q, q') is a two-photon amplitude of the state 
|/3, <?; /3 7 , q 1 ') = bp(q) b^, (g 7 ) |0). For uncorrelated incom¬ 
ing photons Bpp>(q,q') is a direct product of their spec¬ 
tral functions. For two identical photons in the HOM 
geometry, scattering from the TLS with time delay r, 
the amplitude is given by 

Bi2{q,q') = g{q) g(q') e iqT , (n) 


while for the ‘weak blockade’ geometry the appropriate 
function Bn(q,q') differs from the above by symmetris¬ 
ing with respect to q and q'. 

The g-dependent one-photon scattering matrix is given 

by 

~ ( a ’ k \S 1/3, g) = 2 tt S(k - q) s a p(q ), ( 12 ) 

where the on-shell q dependence is equivalent to the w 
dependence given in Eq. 0 . To separate statistical and 
interaction effects, one relates the two-photon S-matrix 
with the appropriate T-matrix, 

s = i-if = i-if°-it int = s 0 -it inU 

where S° (or T°) describes scattering of non-interacting 
photons while *77nt represents the nonlinearity contribu¬ 
tion. The matrix elements of S° are given by 


, k- ak'\S° |/3, g; /3 7 , q') = S%S% + S*S*% (13) 


J ct/3 f u a.'(3'> 


while those of T are found by accounting for the phase 
space restriction on the TLS occupation by photons 


g(ui - w 0 ) = 


/ 1 

(w - w 0 ) 2 

/ /— GXp 

J (JyJ'K 

2 a 2 


but in general all we need to know is its centre, u>q, and 
width, a. As we are considering on-shell scattering, w = q 
(in the units where the group velocity of light is 1 ), we 
may also refer to this function as g(q — wo). 

Then the incoming one-photon state is given by 

l/ 3 )in = J^5(9-w o )&|3(g)|0) . (9) 

When the spectral function centres at |w 0 | = 7 , it cor¬ 
responds to the scattering emulating the balanced beam 
splitter in the HOM geometry, while the zero detuning, 
wo = 0 , corresponds to the exact resonance (or antires¬ 
onance) in the geometry of Fig. [2j In the time-resolved 
representation the photon amplitude, which is a Fourier 
transform of g(q — w 0 ), is a running wave, g{x — t ) cen¬ 
tred at x = t of width of order 1/a. 


(a, k; a', k'\T |/3, g; /3 7 , q') = ^ 2n6(E k -£ q ), (14) 

A-k-A-q 


where 


e k = k + k' , 2 £ k = k-k', = y + * 7 , 

A k = (,k + ij) (.k' + * 7 ) = fi l - i 2 k . (15) 


On substituting the incoming state of Eq. (10), one re¬ 
duces the matrix element Eq. (14) to 

<“'*■" k '\ f \^') =-^r/t $ T ;- (16) 

The probability of finding outgoing photons in the chan¬ 
nels a and a' for the incoming state of Eq. (10 1 is 

f d k d k! 


P n 


1 


'1/3/3' 


N, 3/3' J ( 2 tt )2 


[a,k-,a',k'\ S |/3, /3 7 ) 


3V/3/9' = ( 1 + 5/9/3') J ( 9^2 \ B 0/3’(q^q')\ 2 ■ ( 17 ) 
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Substituting the matrix elements of S, corresponding to 
the independent scattering, Eq. (13), and the nonlin¬ 
earity induced scattering, Eq. (14), we can decompose 
the probability in each channel as the sum P° + SP , 
where SP is the nonlinearity induced change in the prob¬ 
ability of photon bunching (a = a') or anti-bunching 
(a 7 ^ a'). Both in the HOM and resonance geome¬ 
tries we are interested in the anti-bunching probabilities, 
-PhOM = Pl2|12 + P 21 I 12 an d Pres = Pl2|/3^ + -f21|/3/3- 
We focus on the limit T = 7 /cr 1 that corresponds to 
almost monochromatic photons strongly scattered by the 
TLS at the resonance m- In this case the two-photon 


spectral function of Eq. (11) is sharply peaked at q = q' = 
wq while the matrix elements of S' 0 , Eq. (13), are smooth 


functions of wave numbers so that their dispersion can be 
neglected. All further analytic results are found in this 
approximation. 

In the HOM geometry, /? 7 ^ /3', the noninteracting part 
of the probability matrix is governed by the transmission 
and reflection probabilities, T = |f Wo | 2 and R = |r Wo | 2 
with the scattering amplitudes given by Eq. and by 
the following coherence factor, 


v = 


Bi 2 (£k,^k)Bi 2 (£k, —60 , (18) 


whic h is easily found in the Gaussian case, using Eq. 
and (15) as v = e A ~. In these terms we find Phom = 
(T 2 + R 2 ) — 2 vTR. For two identical photons in the 
absence of the delay between them (A = 0) the coherence 
factor v = 1 so that Pjjom = {T — R ) 2 , which vanishes 
when the central photon frequency u/o = 7 (the point 
where the model emulates the balanced beam splitter). 

The nonlinearity-induced correction to Prom is found, 
using Eq. (14) and (Il7|) , as 


<5Phom = —2^12 T 2 (1 — 4R) 


(19) 


where 


C/3/3' - N 00 \ J ^ 


P/3/3'(£fc,Cfe) 


( 20 ) 


Calculating this integral for the Gaussian spectral func¬ 
tion results in £ 12 given by Eq. |?|. 

In the resonance geometry, where both photons are in¬ 
coming from the same channel (j3 = /3 ' = 1 ), one finds 
P r ° s = 2 TR provided that E 1. In this limit in the ge¬ 
ometry of Fig. E h ) independent photons are both trans¬ 
mitted at the exact resonance, T — 1 and R = 1 — T = 0. 
The nonlinearity correction is found similar to that for 
the HOM geometry as 


SP les = Cn T 2 (1 - 4P) 


( 21 ) 


where £n is found from Eq. (20) and In the exact 
resonance P r ° s = 0 and the anti-bunching probability, 
describing the weak photon blockade, reduces to 5P les = 


£11 = £bi- The expression for ^,1 f° r the Gaussian case is 
given by Eq. ([sj) in Results. 

In the opposite limit, T <C 1, the resonance effects are 
weak since the spectral width a greatly exceeds the res¬ 
onance width 7 . This trivial case is of little interest but 
in Fig. [4] we give numerical results for the weak photon 
blockade for intermediate values of larger T up to T = 1 . 
On the inset in this figure we show the full reflection prob¬ 
ability of one out of two simultaneously arriving photons. 
Its maximum at T ~ 1 is higher than (hi due to a trivial 
off-resonance reflection. Note that this result (i.e. in the 
absence of the delay) is in agreement with those earlier 
obtained by Zheng et 


SUMMARY 


In this work we have analysed the role of nonlinearity 
in two-photon propagation in Hong-Ou-Mandel geometry 
and in the resonance geometry in photonic crystals. In 
the ideal HOM geometry two identical photons simulta¬ 
neously arriving at a balanced beam splitter through dif¬ 
ferent incoming ports are bunched together to leave the 
system through the same outgoing port. The HOM inter¬ 
ferometer is thus used for measuring the degree of entan¬ 
glement in biphotons produced by SPDC. Here we argue 
that the presence of the nonlinearity in the beam splitter 
results in partial anti-bunching of the photons thus po¬ 
tentially obscuring such a measurement. We model the 
nonlinearity by a two-level system (which at the same 
time emulates a balanced beam splitter) and calculate 
the anti-bunching probability induced by this nonlinear¬ 


ity. Our results are given by Eq. (19) and (20) and illus¬ 
trated in Fig. [3] 

We use the same model to illustrate the “weak photon 
blockade” of resonant two-photon propagation through a 
TLS. The blockade is due to the impossibility of double¬ 
occupancy of the TLS, and it is always weak due to in¬ 
evitable spontaneous emission. However, exactly at the 
resonance this leads to a finite anti-bunching probability, 
as illustrated in Fig. [4] while two independent photons 
would be either both resonantly reflected (Fig. §0 or 
both resonantly transmitted (Fig. |2jn) . We expect this 
effect to be considerably more pronounced for multiple 
photon propagation which requires further analysis. 
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